Hanbury Brown and Twiss Correlations of Anderson Localized waves 
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When light waves propagate through disordered photonic lattices, they can eventually become 
localized due to multiple scattering effects. Here we show experimentally that while the evolution 
and localization of the photon density distribution is similar in the two cases of diagonal and off- 
diagonal disorder, the density-density correlation carries a distinct signature of the type of disorder. 
We show that these differences reflect a symmetry in the spectrum and eigenmodes that exists in 
off-diagonally disordered lattices but is absent in lattices with diagonal disorder. 



The propagation of quantum-mechanical waves in pe- 
riodic and disordered media is a fundamental theme in 
solid state physics, underlying the transport properties of 
condense matter systems. In a perfectly periodic system, 
the translational invariance gives rise to extended eigen- 
modes known as the Bloch modes. As a result, in peri- 
odic systems an initially narrow wavepacket will expand 
indefinitely and ballistically, i.e. its width will grow lin- 
early in time. Disorder in an otherwise perfectly periodic 
lattice breaks the translational symmetry and can lead 
to exponential localization of the system's eigenmodes 
and to the arrest wavepacket expansion (or diffusion) - a 
phenomena known as Anderson localization 0, ■ 

Traditionally, the localization of waves inside the 
medium was not observed directly, but rather inferred 
indirectly from transmission or conductance measure- 
ments. Recently, a new approach to localization of light 
was realized using disordered photonic lattices 0-0|> m 
which light propagates freely along one axis, and exhibits 
localization in the transverse directions ("transverse" lo- 
calization Q). The equations describing the propagation 
of light in these systems are identical to the equations 
describing the evolution of a single quantum particle in 
an atomic lattice, under the tight-binding approxima- 
tion, thus allowing for the direct observation of Anderson 
localization as originally described in In these ex- 
periments, a localized wavepacket, typically a single site 
wide, was released inside the disordered lattice and al- 
lowed to expand. In periodic lattices, such experiments 
led to the observation of ballistic wavepacket expansion 
0, [l(| ■ As a result of the disorder the wavepacket exhib- 
ited a modified expansion profile, the features of which 
depend on the dimensionality of the system 0, Q , even- 
tually settling to an exponentially localized distribution 
- the hallmark of Anderson localization. Recently, a sim- 
ilar approach enabled the direct observation of Anderson 
localization of matter wave TH- in disordered optical 
potentials, also described by the same equations |14j . 

While it is well established that in Anderson localiza- 
tion the average density distribution exhibits exponential 
localization, not much is known about higher correlations 
of the localized fields. Here we show that spatial intensity 
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Figure 1: (color online). Observation of Anderson localiza- 
tion in lattices with diagonal and off-diagonal disorder. Top 
panels show simulation of the average wavepacket expansion 
in lattices with diagonal (a) and off-diagonal (b) disorder, 
when excited at a single site, showing practically the same 
dynamics, (c), (d) Experimental measurements of the output 
distributions for the two types of disorder (blue lines) as com- 
pared to the same distribution in a periodic lattice (green). 



correlations, also known as Hanbury Brown and Twiss 
correlations [l5| , do not necessarily display a uniform de- 
cay. Moreover, in contrast to the density distribution, we 
find that these correlations carry a signature of the type 
of disorder: while light localized in lattices with on-site 
(diagonal) disorder show decaying correlations, the cor- 
relations in lattices with random tunneling amplitudes 
(off-diagonal disorder) exhibit oscillations. We relate the 
observed correlation features to a spectral symmetry that 
exists in lattices with off-diagonal disorder and is absent 
in lattices with diagonal disorder. Specifically, the eigen- 
values in these disordered lattices are anti-symmetrically 
distributed about the mean value, and eigenmodes as- 
sociated with symmetric eigenvalues share several prop- 
erties. These results are related to a recent prediction 
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of quantum correlations in the single particle limit [l(| , 
and we now show that some analogous features can be 
also be observed in the classical regime. While we shall 
discuss here only the one dimensional problem, these ef- 
fects extend to the two and three dimensional problems 
as well. 

The description for Anderson localization of light in 
one dimensional waveguide lattices is given by set of cou- 
pled discrete Schrodinger equations: 



(a) Diagonal disorder 
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Here n = 1, N where N is the number of lattice sites 
(waveguides), U n is the wave amplitude at site n, fi n 
is the eigenvalue (propagation constant) associated with 
the nth site , C n ^ n ±\ are the tunneling amplitudes be- 
tween two adjacent sites, and z is the longitudinal space 
coordinate (for a more detailed description see e.g. 
These equations are identical to the equations describing 
the time evolution of a single electron in a lattice under 
the tight binding approximation [l| [lfl | , where z repre- 
sent time, and U n is the wavefunction at site n. There- 
fore, while the experiments described in this Letter were 
conducted in the optical domain, the results hold also 
for other systems described by the tight-binding model, 
such as electron in crystalline structure, or Bose-Einstein 
condensate in disordered optical potentials. 

In the tight-binding model, the disorder type falls into 
two broad categories: diagonal disorder, in which the 
(3 n parameters are randomized, but the tunneling ampli- 
tudes C ntn ±i are fixed across the lattice. Such disorder 
was considered by Anderson in his original work in 
what is now known as the Anderson model. With few ex- 
ceptions @, 0| , all the recent experiments reporting the 
observation of Anderson localization of light 0, [B| and 
matter waves [ll, 12 1 were conducted using this type of 
disorder. A second type of theoretically well studied dis- 
order is known as 'off-diagonal' disorder, in which the 
the /3 n parameters are fixed across the lattice, yet the 
tunneling amplitudes are randomized. Such lattices are 
known to exhibit several unique spectral properties 
2l|. However, very little has been achieved so far to ex- 
perimentally observe a signature of these properties. 

The recent experiments on AL of light [3|-|5j and mat- 
ter waves [ill . [l2| , have reported the direct measurement 
of the main features of localization, namely the cross- 
over from ballistic transport to localization as a function 
of time and the level of disorder. For example, Fig. 1 
shows simulations of the wavepacket dynamics in disor- 
dered one dimensional lattices This dynamics starts 
with a ballistic expansion of the wavepacket, similarly 
to the expansion in perfectly periodic lattices. After 
some propagation, a localized component emerges near 
theorigin, co-existing with the transient, ballistic compo- 
nent. As the waves propagate, the ballistic component 
decays and the intensity distribution becomes exponen- 



ts 




-5 5 
(c) Off-diagonal disorder 



cr 



S 




-4 -2 2 4 
Distance 



Figure 2: (color online) Measured density-density correlations 
r r , 9 =< Irlq > I < Ir >< I q > for localized wavepackets in: 
(a) lattices with diagonal disorder. The strong diagonal fea- 
ture reflects finite coherence length of the waves (see text), (b) 
The extracted correlation function, (c) Density correlations in 
lattices with off-diagonal disorder, showing checker-like cor- 
relations, (d) The extracted correlation function, showing 
oscillating correlations. 



tially localized in space. Fig. 1 compares this evolution in 
lattices with diagonal (a) and off-diagonal (b) disorders, 
showing practically identical evolution. 

We have fabricated two types of disordered lattices of 
waveguides, similar to those used to demonstrate Ander- 
son Localization of light In one array, the waveg- 
uide widths were identical, but their separations were 
randomized, modeling off-diagonal disorder. The other, 
a diagonally disordered array, had random- width waveg- 
uides yet constant separations. Light was launched into 
individual waveguides and the output distribution was 
recorded. When the output intensity was averaged over 
~100 realizations (by launching the light at different in- 
put locations) we reproduced the exponentially decaying 
localization shown in Figure 1 (c) and (d) . It is shown in 
comparison with the ballistic expansion that was mea- 
sured in perfectly periodic arrays. We note that the two 
types of disorder led to very similar localized states. 

A significant difference was observed, however, when 
we measured the density correlations of the output dis- 
tribution. Figure 2 presents experimental results of 
density-correlation measurements in disordered lattices. 
Here, for each realization of disorder, the density auto- 
correlation is measured, and then averaged over many re- 
alizations (again by shifting the input site) . The result is 
then normalized, so that T r q =<I r I q > / </,.>< I q >. 
As can be seen in Fig 2, both types of disorder show 
a distinct diagonal feature in the correlation matrix. 
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The length scale of this feature along the main diago- 
nal (q = r) of the matrix is the localization length. The 
width of the diagonal feature is given by the correlation 
length, which is not represented in the ensemble-averaged 
density distributions. In each single realization, the den- 
sity distribution is not a smooth exponentially decaying 
distribution; it is speckled 0, [B| . The width of the di- 
agonal feature reflects the average speckle size. Since in 
each realization the speckles patter varies, their featured 
are smeared out in the averaged distribution, and the in- 
formation about their width is lost. However, the fact 
that these speckles have a characteristic length scale is 
recorded in the averaged correlation function. 

We find that the correlations function carried addi- 
tional information on the type of disorder in the lattice, 
information that is also lost when one considers the av- 
erage density distributions. A closer look at the corre- 
lation matrix reveals that for lattices with off-diagonal 
disorder (Fig 2c) the correlations tend to form a check- 
ered pattern. This can be better seen when looking 
at the correlation function g(Ar), extracted from the 
correlation matrix T r ^ q by summing over the diagonals, 
g(Ar) = r ? . jr+ Ar- The density correlation decays 
smoothly for lattices with diagonal disorder, yet it ex- 
hibits decaying oscillations for lattices with off-diagonal 
disorder. These results were corroborated in numerical 
simulations (not shown). 

To explain the dependence of the density correlations 
on the type of disorder, we start by pointing out a sym- 
metry that exists in periodic lattices that is sustained also 
in lattices with pure off-dia gona l disorder, but not in lat- 
tices with diagonal disorder num. Without loss of gen- 
erality, we can set the diagonal terms of the Hamiltonian 
to zero, and we denote the randomized tunneling terms 
as C n . Let (ai,...,ajv) be an eigenvector with an eigen- 
value A. Then C\ * a% = A * ai ; C± * a\ + Ci * = A * a2 
; C2 * 0,2 + C3 * a 4 = A * 03, etc. Now it is obvious that 
the vector b n = (— l) n * a n , is also an eigenvector with 
an eigenvalue —A. That is, the eigenvectors are paired 
around the center of the band, where twin eigenmodes 
share the same density distribution in absolute value, but 
an opposite (staggered) phase structure. It is also easy 
to see that this property will not be exhibited by lattices 
with diagonal disorder. For a more detailed proof, see 

m 

To visualize this symmetry, we show in figure 3 the 
eigenvalues and eigenmodes for a single realization of a 
disordered lattice with diagonal vs. a lattice with off- 
diagonal disorder. In fig. 3a and 3b, we compare the 
spectrum ("band") of eigenvalues for lattices with N = 50 
sites. As was shown for lattices with diagonal disorder 
0, disordered lattices support two types of tightly lo- 
calized eigenmodes with eigenvalues at the edges of the 
spectral band: at one edge the eigenmodes are tightly lo- 
calized in space, each mode occupies a different location, 
and they generally have a flat phase profile (zero phase 
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Figure 3: (color online). Spectral symmetry in lattices with 
diagonal and off-diagonal disorder, (a) The spectrum (band) 
of eigenvalues for lattices with off-diagonal disorder. Note 
the perfect symmetry with respect to the band center, (b) 
The spectrum for lattices with diagonal disorder, showing no 
such symmetry, (c) Pairs of eigenmodes taken from symmet- 
ric eigenvalues occupy the same region of the lattice, have the 
same spatial distribution, but vary in phase structure, (d) In 
lattices with diagonal disorder no such symmetry exists, (e) 
Experimental measurement of two spectral-symmetric eigen- 
modes in a lattice with off-diagonal disorder, (f) A similar 
measurement in a lattice with diagonal disorder cannot ex- 
cite two modes at the same position (see text). 



difference between adjacent sites). At the other edge of 
the spectrum the modes are also tightly localized, only 
that they are staggered - there is a ir phase difference be- 
tween adjacent sites. The eigenmodes at the center of the 
band have a more complicated phase structure, and they 
are typically wider. As we concluded above, in lattices 
with off-diagonal disorder we find that each eigenmode 
at one edge of the spectrum had a twin-eigenmode at 
the other edge. These twin eigenmodes occupy the same 
region of the lattice, have the same distribution of the 
absolute amplitude, but they differ in phase structure as 
shown in Fig. lc. This property does not exist in lattices 
with diagonal disorder (see Fig. Id). 

This special symmetry of the can be observed exper- 
imentally. Briefly, it is possible to excite pure localized 
eigenmodes by using either a flat-phase or a staggered 
beam, with the correct width and initial position at the 
lattice. Flat-phased and staggered localized eigenmodes 
of a lattice with diagonal disorder were measured (See 
figure 2 in Q), showing indeed that in diagonally disor- 
dered lattices the positions of the lowest-eigenvalues flat- 
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phased eigenmodes never coincided with the position of 
the highest-eigenvalues, staggered localized modes. Fig. 
3f shows an excitation of a flat-phased eigenmode in a lat- 
tice with diagonal disorder. Indeed, when the same input 
beam was tilted to excite neighboring sites with a ir phase 
difference, the output density showed considerable expan- 
sion, suggesting that no staggered localized eigenmode 
resides in the same location. In contrast, Fig 3e, shows 
the same procedure in a lattice with off-diagonal disor- 
der. Here, the flat phased beam excited a flat-phased, 
localized eigenmode, and a beam with a ir phase differ- 
ence between adjacent sites excited a staggered localized 
eigenmodes with the same width, the same spatial profile, 
and at the same location. 

Now, to explain the different density-density correla- 
tions shown in figure 2, we need to consider the effect 
of these different spectral properties on the expanding 
wavepackets when they are excited at a single site. The 
localization effect arises from the fact that all the eigen- 
modes that are excited have a finite extent. It is well 
known that in infinite one dimensional disordered sys- 
tems all eigenmodes of the system are localized. In lat- 
tices with off-diagonal disorder, an initial excitation of 
a single lattice site necessarily involves the simultane- 
ous excitations of pairs of "twin" eigenmodes, as they 
have identical overlap with the initially excited site. The 
sum of two identical amplitude distributions yet with a 
7r-phase difference in each second site results in a den- 
sity comb-like pattern that nulls at every second site. 
In the dynamic problem several pairs could be excited 
simultaneously by the single site initial condition, and 
the two modes of the pair accumulate phase in a differ- 
ent rate (according to their eigenvalues). Nevertheless, 
the wavepacket will contain a component with an oscil- 
lating intensity pattern, with a spatial frequency of two- 
sites. This effect is washed out in the density distribution 
averaged over all realizations of disorder, as in each re- 
alization the oscillations appear in a different location. 
However, the fact that such oscillations appear in each 
realization will be recorded in the averaged correlation. 

We note a previous publication by our group [16], 
in which oscillating Qitanfumcorrelations were analyzed 
theoretically for bosonic or fermionic pairs, predicting 
checker-like correlations in some cases. Those phenom- 
ena are not unrelated to the results reported here, yet 
the oscillating correlations reported there were of quan- 
tum origin (i.e. in the case of light they require the use of 
non-classical light), while here the described effects are 
purely classical, wave effects. 

In conclusion, we have experimentally shown that den- 
sity correlation measurements can carry a signature of 
the type of disorder that exists in a given sample, and we 
have traced that signature to the existence of a unique 



spectral symmetry that is exhibited by lattices with off- 
diagonal disorder. Similar results can be measured in 
matter waves systems, using density-correlation measure- 
ments [24] in disordered lattices. It might also be in- 
teresting to study the effect of nonlinear interactions on 
these correlations, either in the optical or matter wave 
experiments. The signature of these results might also 
be observed in correlations measurements for multiply 
scattered classical [26j and non-classical [27, 28 1 light. 
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